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I. INTRODUCTION 

Nowadays the standard model (SM) becomes trustworthy to describe microscopic funda¬ 
mental physics, since the SM Higgs has been discovered at the CERN Large Hadron Collider 
(LHC). However it has to be still extended in order to include dark matter (DM) and (mi¬ 
nuscule but) massive neutrinos that make an allusion to their existence by overwhelming 
experimental evidences. One of the economical and elegant solutions to resolve these con¬ 
troversial issues is to accommodate radiative seesaw models to SM, in which active neutrino 
masses are radiatively arisen and exotic helds are naturally introduced to induce such tiny 
masses. Such an exotic held can be frequently identihed as a DM candidate. In this sense, 
one might hnd that neutrinos have a strong correlation to the DM candidate. Following 
the landmarks HHZI. a vast of literature on radiative seesaw model has recently arisen in 
Refs. [8H92] 

In this paper we employ new fermions and bosons in addition to the SM-like Higgs boson, 
in which the leading neutrino masses can be induced at the two-loop level where the relevant 
Lagrangian is controlled by an additional global U{1) symmetry. And the effective tri-linear 
coupling between the SM-like Higgs and an isospin triplet boson, which is needed to have a 
massive CP-odd neutral boson, is also generated at the one-loop level through such exotic 
helds after the global U{1) symmetry spontaneously. 

Following the paper of [^, the neutral component can be a DM candidate if there is 
enough mass diherence between two neutral fermions to evade the constraint of the direct 
detection searches via the SM neutral gauge boson (Z). The mass diherence arises from the 
type-H like term after acquiring the vacuum expectation values (VEVs) of the isospin triplet 
boson. However since the mass diherence is very tiny because its VEV be less than a few 
GeV constrained by the electroweak precision test, we have to work on the co-annihilation 
system to obtain the observed relic density, if we focus on rather lighter DM mass that is 
less than 80 GeV. 

Heavy charged lepton masses themselves are constrained by the LEP and LHC experi¬ 
ment, and the mass diherence between the DM mass and its lepton is also constrained by 
the electroweak precision test. As a result, the allowed range of the DM mass can be highly 
restricted. 

We have two new sources to explain the deviation of anomalous magnetic moment to 
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SM. However either of them that has a strong correlation to the neutrino masses cannot 
reach the sizable value of the anomalous magnetic moment. This is because the neutrino 
oscillation requires rather large off-diagonal neutrino mass matrix elements, which tends to 
be in conflict with the anomalous magnetic moment. Notice here that constraints of the 
lepton flavor violations (LFVs) that always emerge in such radiative models are not so strong 
to restrict the Yukawa couplings related to the neutrino masses. 

This paper is organized as follows. In Sec. II, we show our model building including Higgs 
masses, neutrino mass, LFV, muon anomalous magnetic moment, and DM. In Sec. HI, we 
show our numerical results. We conclude in Sec. VI. 

II. THE MODEL 


Fermion 

Ll 

efi 

T f 

^L{R) 

^'l(r) 

{SU{2)l,U{1)y) 

(2,-1/2) 

(1,-1) 

(2,-1/2) 

(1,-1) 

0(1) 

-1 

-1 

-3/2 

-3/2 


TABLE I: Lepton sector; notice the three (or two) flavor index of each field L^, br, ^'l{r) 
®L(R) abbreviated. 


Boson 

$ 

V 


A 

S 

(5C7(2)i.,[/(l)y) 

(2,1/2) 

(2,1/2) 

(1,0) 

(3,1) 

(1,0) 

0(1) 

0 

3/2 

1 

3 

1/2 


TABLE 11: Boson sector 


In this section, we review our model, in which the particle contents for leptons and bosons 
are respectively shown in Tab. 
isospin doublet and with isospin singlet, but —3/2 charge under the global symmetry 
to the SM fields. Each of the exotic field needs (at least) two flavors in order to satisfy 
current neutrino oscillation data [9l]. Moreover, an introduction of e' is requested only to 
have mass difference between the neutral component and the charged component of L' in 
our model. Otherwise our DM (neutral component of L') does not satisfy the observed 
relic density VLh? ?«0.12 [95] due to an enhancement of its thermal averaged cross section 


and Tab. 


0 


We add vector-like fermions of with 
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with co-annihilation system, as we will discuss later. As for new bosons, we introduce two 
neutral isospin singlet scalars ip and S with 1 and 1/2 global charge for each, an isospin 
doublet scalar rj with the 3/2 global charge, and an isospin triplet scalar A with the 3 global 
charge. Note here that <h is neutral under the global charge not to couple to our physical 
Goldstone boson. Then we assume that <h, 99 , and A have VEV, which are symbolized by 
n/v^, v'/\/2, and va/ respectively. VEV of ip spontaneously breaks the global symmetry 
down. Even after the global U{1) breaking as well as the electroweak breaking, a remnant 
discrete symmetry Z 2 remains, which is understood as an accidental symmetry. This Z 2 
symmetry plays a role in assuring the stability of our DM candidate; neutral component of 

j' 

^L(Ry 

The relevant Lagrangian for Yukawa sector and scalar potential under these assignments 
are given by 

-Cy = {ye)aLL,^eR^ + {yL)aiLL,L'^^S + + {y'j^)ijL'^.{iT2)AL'j^^ 

+ iy'LR)ijL'L,^e'^j + iy'RL)ijL'R^^eL. + {ys)ibe'L^eR^S* 

+ + h.c. (II.l) 

V = -h rri^\ri\‘^ + + m^Tr[|Ap] 

+ h.c.) -h {iT 2 )A^r] + h.c.) -h Xoiy'^^Sip + h.c.) 

+ + A'^J$t77|2 + A<i,^|$| Vr + VP + A<i,A|$pTr[|Ap] 

1-3 

+ ^(<hVj<h)Tr[AViA] -h A,,|?7p -h A,,^|?7| VP + VI^Pl'^'P + Ar,A|hpTr[|Ap] 

i 

+ V 5^(WTr[AViA] + Wp>\^ + A^5|(^p|^p + A<^A|<^pTr[|Ap] 

i 

+ A 5 |^P + A 5 A|^pTr[|Ap] + AA(Tr[|Ap])2 + X'^Det[A^A], (II.2) 

where Ti{i=l-3) is Pauli matrix, each of the index a{b) and i{j) that runs 1-3 and 1-2(3) 
represents the number of generations, and the hrst term of Cy can generates the (diagonal¬ 
ized) charged-lepton masses. We work on the basis where all the coefficients are real and 
positive for our brevity. 
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A. Scalar sector 


After the EW symmetry breaking, each of scalar held has nonzero mass. We parametrize 
these scalar helds as 


$ = 


kO 


, V = 


7 ]'' 

rj 


,0 


A = 


^ A"" 
-7i 


(II.3) 


And the neutral components of the above helds and the singlet scalar held can be expressed 
as 

0° = + h + m), rf^ = + ir//), A° = + Ar + iA/), 


V= 4(”'+s=4(s«+ ^Si), 


where h is the SM-like Higgs, and v and ua is related to the Fermi constant Gp by v'^ + 2v\ = 
1/(v^Gr) (246 GeV)2. 

The CP even Higgs boson mass matrix with VEV in the basis of (Ar, h, p) is given by 


r2\vev 


(M ;Qp_gygn 


/ 


V 


+ 2AA'yA [(A$A + A^a)'^A ~ V^Aieff] V XAifiV'vA 

[(A$a + A$a)'^a - V2pefi] V 2A$n^ 


\ 


Xa^v'va 




2X^v 


/ 
(11.5) 


and (M^)((fp_gygj^ is diagonalized by 3 x 3 orthogonal mixing matrix Or as 
OR(M2)^f_g^gj^OR = diag.(m^^,m^gj^,m^ 3 ). Here hsu is the SM Higgs and hi and hs are 
additional Higgses, eigenstate. 

The CP odd Higgs boson mass matrix with VEV in the basis of (A/, a) is given by 

~'^heSV \ 

2v^P'efr^A J 

and {M‘^)cp-odd i® diagonalized by 2 x 2 orthogonal mixing matrix Oj as = 

diag.(0, ml), where ml = massless mode is absorbed by the neutral gauge 

boson Z to be massive. 

The CP even inert Higgs boson mass matrix in the basis of (?7r, S'r) is given by 


(M' 


2 \i;ef 

CP-odd 


( 11 , 6 ) 


(M 


2\inert 
/CP—even 


2m^ 

2 2 
\„vv' 2m|+2P2Atsi;'+Ay,gi;'^+Aj.g-i)^+AgAA 


(11.7) 
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and is diagonalized by 2 x 2 orthogonal mixing matrix Vr as 

VR{M‘^)^fX.eryR = diag.(m^^,m^J. 

The CP odd inert Higgs boson mass matrix in the basis of (r^/, Sj) is given by 


(M^) 


inert 

CP-odd 


2 

Xqvv' 

2 


Xqvv' 

2 

2m|-2p2/xsi)'+A^Si)'2+A^5i;2+AsA''^A 

2 




( 11 . 8 ) 

and is diagonalized by 2 x 2 orthogonal mixing matrix Vr as = 

The singly charged Higgs boson mass matrix with VEV in the basis of is given 


by 


(M 


2\vev 
1 singly 


iV2fi^fi-X'^^VA)v'^ f X'^,^VA _ 

2va \ V2 heSy 


)■ 


(II.9) 


and (Tf^)gj^giy is diagonalized by 2 x 2 nnitary mixing matrix Oc as Oc(M^)g^gjyO^ = 
diag.(0,m^), where ttiq = +2^a) massless mode is absorbed by the 

charged gange boson to be massive. The singly charged inert boson mass eigenstate is 
given by 


m^± = 


2ml + + (A^a - A^a)^z 


( 11 . 10 ) 


The donbly charged boson mass eigenstate is given by 


2 _ ( /^eff 


<S>A 


( 11 . 11 ) 


B. Effective trilinear coupling of //eff 

In onr model, the term /iefr<h^(ir 2 )Al$ is forbidden at the leading order, bnt indnced at 

the one-loop level mediated by inert nentral bosons r]^ and S as depicted in the lower part 

of Fig. [T] The effective trilinear conpling of fies is given by 

(Aon')V»? dxidx2dx3S{xi + X2 + X3 - 1) m 

hefF ~ a I (A \2 

6 Jo m 

\v^)2^{VR)^l{V^)lAVR),2{V^)lk{VR)^ _ (H/)2.(H,)a(H/)i,(V}),2(ff)u(V}),i l ^ 

+ ^2mH. + xim\_ -f X2m\. + X3m\^ 

where each of {i,j, k) rnns form 1 to 2. 
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C. Inert conditions 


To forbid VEVs for our inert bosons and S, the quartic couplings of \ and \s has 
to be always positive. To achieve the situation up to one-loop level, we have to satisfy the 
following conditions at least: 


0 < A“*' = A, + (5A<‘I, 0 < Af *' = As + (5A^">, 


(11,13) 


where 


AA!‘> = -10KI‘‘f„(),,A), 


(11,14) 


AA™ = -8|,^sl‘‘ro(S,^), 


(11,15) 


ro(/i,/2 ) = tjT f dxdyi(x + y-l)- - 

(47r)2 Jo [xmj^ + ymj^ 


( 11 . 16 ) 


D. Fermion Sector 


Let us hst dehne the exotic fermion as follow: 


T' — 
^L(R) - 


(11,17) 


Neutral exotic fermion: Then the mass matrix for the neutral fermion in the basis of [Nf, Np 
is given by 


Mm — 


Ml m 
m Ml 


( 11 . 18 ) 


where m = {m « Ml) and we assume to be one generation case with positive 

real for our simple analysis. M^ is diagonalized by 2 x 2 unitary mixing matrix Viv as 
VnMmV^ = diag.(Mx, — m, Ml + m), where Vn is written as a maximal mixing form 


1 i —i 

V2 1 1 


( 11 . 19 ) 
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and we define the mass eigenstate of the nentral fermion [A^i, iV|] then we have the following 
relations: 


- 1 

1 

iNi + Ntf 

1 


—iNi A- N 2 


( 11 . 20 ) 


Fnrthermore we redefine these helds as 'ijji = Ni + and ^’2 = -^2 + -^ 2 , then we obtain 
the Majorana helds ipi and '02- We snmmarize the relations between -0 and N below 


PiV'i = W, PrA = ivr, Pl^2 = N^, Pr^2 = N 2 , (11.21) 

'iP.Pr = w, = iV{=, tP^PR = m. ^2 Pl = N 2 . (11.22) 


The lighter held Ni can be a DM candidate, bnt the mass diherence between them 2m is 
expected to be tiny becanse m is originated from va- As a result, we have to consider the 
co-annihilation system (at least) including t/ji and 1^2 to obtain the relic density. Notice here 
that the lowest bound on m comes from the inelastic scattering through Z boson and should 
be heavier than 0(100) keV [UH] . 

Singly charged exotic fermion: The mass matrix for the singly charged fermion in the 
basis of [E', e'] is given by 


Me = 


Ml m'g 
m; Mr 


(11.23) 


where we assume m'^ = yfj^v/\/2 y'j^]^v/\/2 (m'g << MrMr) with one generation. Note 

here that the mass matrix can be taken positive real without loss of generality. Mr is 
diagonalized by 2 x 2 unitary mixing matrix Vc as 

VcMeV^ = ^diag. {Ml + Mr- V(Mi-M^)2 + 4m;, Ml + Mr+ ^{Mr - MrY + Am')j , 

(11.24) 


where Vc is as 


w = 


Ce Se 
—se Ce 


m' 


S2E — 


(11.25) 


^{ML-MRy + AmP 

and Se = sin^E and ce = cos 6 e- Then we dehne the mass eigenstate of the charged fermion 
[Pi, P 2 ], we have the following relations: 


E' 


CrEi — SeE2 

e' 


SrEi -\- CeE2 


(11.26) 
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Relation between neutral and charged fermion: Since the DM mass {Mx) can be approx¬ 
imately given by and m'^ « we have the relation from the charged fermions 

approximately: 

Me, ^ Mx, Me, ^ Mr. (11.27) 

On the other hand if the mass difference between X and Ei is enough tiny, Ei also partic¬ 
ipates in the co-annihilation system of DM. The upper bound on Me is derived from the 
electroweak precision data, which is typically written in term of AS, AT, AU parameters. 
The most stringent bound comes from AT that suggests the mass difference between Mx 
and Me should be smaller than 45 GeV [93] : 

Me, - Mx < 45 GeV. (11.28) 

The mass of Ei is also constrained by the LHG searches and the lower bound is around 
300 GeV if Ei decays into the SM charged-lepton -|- missing(=DM), which is an analogous 
analysis of the slepton searches |S7|. However Ei cannot decay into the SM charged-lepton -|- 
missing kinematically but decay into the SM charged lepton -|- two missing fermions( active 
neutrinos -I- DM) through the lighter mass of Hi or Ai, which is assumed to be heavier than 
the mass of Ei in our case. Then only the constraint comes from the LEP experiment that 
suggests [96] 

100 GeV < Me, and 15 GeV < Me, - Mx- (11.29) 


E. Neutrino mass matrix 


The neutrino mass matrix can be generated at two-loop level as depicted in Fig. [TJ and 
its form is given by 

T 1,2 1,2 r 




2(47r) 


'^iyL)a,k{yL)b,k ^ 




T ',2 / 
2i> 


^1,. - 1 


lnXi,j 


X. 


2,i 


InXsj 




T ^2 / 


R,2j, 


.Y 2 i 

1 ^- 

Yi,j - 1 Y2,j - 1 


(11.30) 


where Xjj = and Xij = ' Remind here that two flavor of Efc(k=l- 

2) is introduced to obtain the current neutrino oscillation data. {Aiu)ab can be generally 
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L' 




5 


V 



I'L 


FIG. 1: Feynman diagram for the neutrino mass. 


diagonalized by the Maki-Nakagawa-Sakata mixing matrix Vmns (MNS) as 

C13C12 C13S12 sise-*"^ 


14dNS 


—C23S12 — S23Sl3Ci2e*'^ C23C12 — 523^135126*^ 523C13 5 

523512 — 6235136126®'^ —523C12 — 6235135126®'^ 623613 


(11.31) 

(11.32) 


where we neglect the Majorana phase as well as Dirac phase 6 in the numerical analysis for 
simplicity. The following neutrino oscillation data at 95% conhdence level [M] is given as 


0.2911 < 5^2 < 0.3161, 0.5262 < 5^3 < 0.5485, 0.0223 < 5^3 < 0.0246, (11.33) 

\ml^ -ml^\ = (2.44 ± 0.06) x lO’^ eV^ - ml^ = (7.53 ± 0.18) x 10"^ eV^ 


where we assume one of three neutrino masses is zero with normal ordering in our analysis 
below. 


F. Muon anomalous magnetic moment and Lepton flavor violations 

The muon anomalous magnetic moment (muon g — 2) has been measured at Brookhaven 
National Laboratory. The current average of the experimental results is given by m 


= 11659208.0(6.3) x lO"'”. 
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Process 

ihJ) 

Experimental bounds (90% CL) 

y- e-y 

(2,1) 

Br(/i — ;■ ej) < 5.7 x 10“^^ 

T~ — >■ e“7 

(3,1) 

Br(T —>• ey) < 3.3 x 10“® 

r" ^ /i"7 

(3,2) 

Br(r —)• py) < 4.4 x 10“® 


TABLE III: Summary of li —)> ij'j process and the lower bound of experimental data [10U| . 


It has been well known that there is a discrepancy between the experimental data and the 
prediction in the SM. The difference Aa^ = — a™ was calculated in Ref. [HE] as 


Aa^ = (29.0 ±9.0) x 10“^^ 


(11.34) 


and it was also derived in Ref. |99] as 


Aa^ = (33.5 ±8.2) x 10 


-10 


(11.35) 


The above results given in Eqs. (11.34) and (11.35) correspond to 3.2 ct and 4.la deviations, 
respectively. 

In our model, we have new contributions to Aa^ coming from yi and ys terms. These 
contributions are calculated as 


Attf, ^ 


E Bl) + Bl)] 

^ ' a 

(vaV rr?2 

E [(±)LG(ffa. E,) + E^)] 


where 


G(/i,/2) 


^1 nl—x 

dx I dy 


x{y + z) 


'0 


0 xm\ + {y + z)m\ 


z=l—x—y 

4 ^2 


2m ± Smj-^mj^ — ± m^^ ± In 


(11.36) 


(11.37) 


6(mj^ — 

Our relevant lepton flavor violation process {ii —)■ ijj) comes from the same terms of 
anomalous magnetic moment at the one-loop level in principle. Each of flavor dependent 


process has to satisfy the current upper bound, as can be seen in Table III However the 


contribution from ys can be always negligible assuming the diagonal ys- This is because 
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this term does not contribute to the neutrino masses. Hence we consider the contribution 
from ui only. Then the branching form is given as 


1-2 


Br(£i £j7) cs ^ 


^\{yL)jk{yL)ik\ 


(167r2GF)2 


)ia9{Ha. E,) + (H/ E,) I , (11.38) 


where 

9{f 1 J 2 ) 

and Gp is Fermi constant. 


nl—x 

dx / dy 


xy 


g xm‘j_^ + (?/ + 


/2 


(11.39) 


z=l—x—y 


G. Dark matter 


First of all, we discuss the direct detection searches reported by the experiment of 
LUX nnn. As mentioned before, the inelastic scattering process through Z boson is al¬ 
ways evaded to retain the mass difference with 0(100) keV, which is generated via y'j^VA in 
our model. We have Higgs portal process and our scattering cross section with nucleon is 
given by 


(Ttv ~ 0.082- 


m 


nWl 


71 


1-3 

X 


(oS..(oS« 


mt 


(11.40) 


where h 2 = hsu and the mass of neutron, which is symbolized by mN, is around 939 GeV. 
LUX suggests that should be less than 0(10“'^®) cm^ at 0(10) GeV mass range of DM. 

Relic density of DM: Before the serious analysis to the relic density, we £x a situation on 
DM(which is denoted by -01 = X) as follows. We focus on the region of 50 GeV < Mx ^ 
80 GeV to realize a hne perspective(or simple) analysis, and suppose a co-annihilation system 
on DM, in which there exist degenerated fermions "02 and Ei that could affect to the relic 
density (Dh^ ~ 0.12) reported by Planck [95]. Its complete analysis and formula can be 
found in Ref. HE] and HE]. Our relevant processes for the thermal averaged cross section 
comes from the co-annihilations and annihilation for the degenerated helds of 4/j lUi —)■ £i>l 
via W boson, and EiEi —>■ // via t-channel mediated by Ha and Aa that comes from y^ and 
s-channel mediated by Z boson. Here / represents all the SM fermions that are kinematically 
allowed. As subdominant modes, there exist XX —)• via t- and u-channels mediated by 

Ha and A^, X'02 —)■ // via s-channel mediated by Z boson, XX ^ 2G via t- and u-channels 
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mediated by 1 ^ 2 , and XX —)■ // via s-channel mediated by three CP-even neutral bosons hi 
from y'j^. Once the total cross section is given in al the above processes, the effective 

annihilation cross section is given by 

k k 

^eff^rel EE (1 + A.)®-'yi + (11.41) 

i=l j=l 

where ^^efr is the effective degree of freedom 

k 

= + (11.42) 

i=l 

and Aj = [rrii — mi)/mi is the mass difference between DM and the other degenerate 
particles, gi = 2 is the degree of freedom for each Majorana particle Xi) ^ = ^i/T and aijV 
is (co-)annihilation cross section between i and j. Then the formula of relic density is found 
as 


VLh^ 


1.07 X 10^ 


jl'^AfpilGeVI r (ijf + eta) ’ 


(11.43) 


where ~ 100 is the total number of effective relativistic degrees of freedom at the time of 
freeze-out, Mpi = 1.22 x 10^®[GeV] is Planck mass, x/ ~ 25, and Oefr and Oeff are derived by 
expanding UefrUrei in terms of Urei as 


^effDel ~ ®eff T ^reTeflf" 


(11.44) 


Notice here that the mass of DM is assumed to be less than the mass of boson(^ 81 
GeV), because the cross section is too large to satisfy the relic density once the hP^-boson 
hnal modes are open. 


III. NUMERICAL RESULTS 

Now that all the formulae have been provided, we have a numerical analysis. First of all, 
we £x the following parameters in the scalar sector: 

UA ~ 0.0007686 [GeV], v' ^ 319.7 [GeV], ys ~ 278.0 [GeV], ~ 715.2 [GeV], m^ ^ 371.8 [GeV], 
ms ^ 341.6 [GeV], A$ 0.1300, Aai ^ 0.7226, Aa^ ^ 0.3587, A,, ^ 0.9717, A^ ^ 0.6274, 

As ^ 0.9333, A$a ~ 0.8641, A'^^ = 0.7313, A,,a - 0.9876, ^ 0.7700, A^a ~ 0.2470, 

Asa ~ 0.7253, ^ 0.6639, A'^,, ^ 0.9232, A<i,^ ^ 6.912 x 10"®, As$ ~ 0.7578, 

A^^ ^ 0.3144, As^ ~ 0.9915, As<^ ~ 0.6321, Aq ^ 0.6693, /ies ^ 0.002806 [GeV], (III.l) 
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where these above values satisfy the inert conditions in Eq. (11.13). Then we obtain the 
physical values as follows: 

m/ig^ ~ 125.5 [GeV], rrih^ cs 368.1 [GeV], irihi ~ 395.3 [GeV], 

iriHi ~ 442.2 [GeV], ~ 551.7 [GeV], niAi ~ 204.9 [GeV], tua^ ~ 454.6 [GeV], 


rUr 


O 


R 


-1 7.4 X 10-® 2.1 X 10-^ 


r ~i 

-7.4 X 10“® -1 4.8 X 10“® 

, Oj = 

6.2 X 10“'’ 

1 

-1 

6.2 X 10-6 

2.1 X 10-^ 4.8 X 10-® 1 




Vr 


-0.97 0.25 
0.25 0.97 


, V, 


-0.16 0.99 
0.99 0.16 


, Gc 


4.4 X 10-® 1 

-1 4.4 X 10“® 


. (III.2) 


We search the other physical values with these values, where we take 

0.9 <yL<V2, 10“® < ( 2 /l) 2,3 < 10■^ 50 GeV <Mx <80 GeV, 
Mei < Mx + 45 GeV, 300 GeV < *,] ^ 500 GeV, 


(III,3) 


in our numerical parameter spaces to reproduce neutrino oscillation data, LFVs, and relic 
density of DM by using these above values. Notice here that we solve the other 


by comparing the theoretical form and the central experimental value in Eqs. (11.30)-(11.34), 


and Me 2 is assumed to be large enough to decouple these phenomenologies . Here we show 


four benchmark points (BPs) as shown in Table IV, in which we obtain Br(r —)■ ey) and 
Br(r —)■ fi'y) for all the BPs are respectively 0(10“^^ ~ 10“^^) and 0(10“^^ ~ 10“^^) that 
are completely safe at the current bounds. On the other hand Br(/i —ey) is rather close 
to the current upper bound. However the BPi and BP 2 does not satisfy the LEP bound in 


Eq. (11.29), since the DM mass is less than 100 GeV while the mass difference between Mx 


and Mei is greater than 15 GeV. As a result, only the BP 3 and BP 4 are complete solutions 
to satisfy all the data that we discuss. 

As for the direct detection searches, our elastic scattering cross section form reduces to 
1.1 X 10~^^{y'i)‘^ [cm^] at the range of O{10) GeV DM mass, which is always below the 
current upper bound of the LUX experiment. 

As for the anomalous magnetic moment, Aa^ = 0(10“^^) at most is obtained by the 


contribution from yi in the hrst term of Eq. (11.38), which is much below the sizable value 
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l(yL)i,il 

(yi)i ,2 

(yi)i,3 

\{Vl)2,i\ 

( 2 /^) 2,2 

(yi)2,3 

Mx 

GeV 

Me^ 

GeV 


Br(/i —> eq) 

BPi 

0.0025 

0.031 

0.036 

1.6 X 10"^ 

2.7 X 10-^ 

4.5 X 10-5 

62.27 

79.55 

0.119 

8.1 X 10-^^ 

BP 2 

0.0029 

0.028 

0.032 

2.5 X 10-5 

3.9 X 10-5 

4.3 X 10-5 

59.65 

76.49 

0.120 

9.6 X 10-^^ 

BP 3 

0.0023 

0.029 

0.034 

5.4 X 10-5 

9.5 X 10-5 

4.5 X 10-5 

78.67 

103.6 

0.112 

7.3 X 10-^^ 

BP 4 

0.0014 

0.0216 

0.019 

2.9 X 10-5 

5.1 X 10-5 

4.7 X 10-5 

79.7 

100.5 

0.125 

7.8 X 10-^5 


TABLE IV: Four bench mark points(BPs), where Br(r —)• 67 ) and Br(r —)• ^'y) for all the BPs are 
respectively 0 ( 10 “^^ ~ 10 “^^) and 0 ( 10 “^^ ~ 10 “^^) that are completely safe. 


in Eq. (11.34) or Eq. (11.35) . It originates from the fact that yi cannot enlarge due to the 
the constraints of neutrino oscillation data and LEV of /r —)■ eq. On the other hand the 


contribution from ys in the second term of Eq. (11.38) reaches the lower experimental value 
in the perturbative limit(i.e., ys = dvr), since this term can be independent of such kind of 
constraints. In this limit we obtain, for example, Aa^ k. 2 x 10“® with Me 2 =600 GeV. 


IV. CONCLUSIONS 

We have studied a two-loop induced radiative neutrino model, in which we have shown 
some allowed bench mark points to satisfy the observed neutrino masses, LFVs, and the 
relic density of DM in the co-annihilation system satisfying the current upper bound on 
the spin independent scattering with nucleon as well as LEP. We have also shown a new 
source {ys) to marginally obtain the sizable value of the anomalous magnetic moment in the 
perturbative limit. 
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